Draft version January 19, 2010 

Preprint typeset using LAT^X style cmulatcapj v. 04/20/08 



DISCOVERY OF A NEW SOFT GAMMA REPEATER: SGR J0418 + 5729 

A. J. VAN DER HORST 1 , V. CONNAUGHTON 2 , C. KOUVELIOTOU 3 , E. GOGU§ 4 , Y. KANEKO 4 , S. WACHTER 5 , M. S. BRIGGS 2 , 

J. Granot 6 , E. Ramirez-Ruiz 7 , P. M. Woods 8 , R. L. Aptekar 9 , S. D. Barthelmy 10 , J. R. Cummings 10 , 
M. H. Finger 11 , D. D. Frederiks 9 , N. Gehrels 10 , C. R. Gelino 12 , D. M. Gelino 13 , S. Golenetskii 9 , K. Hurley 14 , 
H. A. Krimm 1 ", E. P. Mazets 9 , J. E. McEnery 10 , C. A. Meegan 11 , P. P. Oleynik 9 , D. M. Palmer 15 , V. D. Pal'shin 9 , 
A. Pe'er 16 , D. Svinkin 9 , M. V. Ulanov 9 , M. van der Klis 17 , A. von Kienlin 18 , A. L. Watts 17 C. A. Wilson-Hodge 3 

Draft version January 19, 2010 

ABSTRACT 

On 2009 June 5, the Gamma-ray Burst Monitor (GBM) onboard the Fermi Gamma-ray Space Tele- 
scope triggered on two short, and relatively dim bursts with spectral properties similar to Soft Gamma 
Repeater (SGR) bursts. Independent localizations of the bursts by triangulation with the Konus-RF 
and with the Swift satellite, confirmed their origin from the same, previously unknown, source. The 
subsequent discovery of X-ray pulsations with the Rossi X-ray Timing Explorer (RXTE), confirmed 
the magnetar nature of the new source, SGRJ0418 + 5729. We describe here the Fermi/ GBM ob- 
servations, the discovery and the localization of this new SGR, and our infrared and Chandra X- 
ray observations. We also present a detailed temporal and spectral study of the two GBM bursts. 
SGR J0418 + 5729 is the second source discovered in the same region of the sky in the last year, 
the other one being SGR J0501 + 4516. Both sources lie in the direction of the galactic anti-center 
and presumably at the nearby distance of ~ 2 kpc (assuming they reside in the Perseus arm of our 
galaxy). The near-threshold GBM detection of bursts from SGR J0418 + 5729 suggests that there 
may be more such "dim" SGRs throughout our galaxy, possibly exceeding the population of "bright" 
SGRs. Finally, using sample statistics, we conclude that the implications of the new SGR discovery 
on the number of observable active magnetars in our galaxy at any given time is < 10, in agreement 
with our earlier estimates. 

Subject headings: pulsars: individual (SGR J0418 + 5729) — stars: neutron — X-rays: bursts 



1. INTRODUCTION 
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In the last decade, observational evidence for neutron 
stars with extreme surface dipole magnetic fields (B ~ 
10 14 - 10 15 G) or "magnetars" ()Duncan k Thompson! 
1992) has steadily grown, with more than 15 magnetar 
candidates to date. The majority of magnetars are mem- 
bers of two neutron star populations historically known 
as Soft Gamma Repeaters (SGRs) and Anomalous X-ray 
Pulsars (AXPs); a couple were previously classified as 
Isolated Neutron Stars or Compact Central Objects. Al- 
though these systems have tangible differences, they are 
also linked with a multitude of similar properties (for re- 
cent reviews, see iWoods fc Thompsonll2006l: iMereghettH 
[20081) . Most sources are visible only in X- and low-energy 
gamma rays; very few have also been detected in the op- 
tical and infra-red, w hile two sources have been observed 
at radio wavelengths (jCamilo et al.ll2006L 120071 ) . All but 
two reside in our Milky Way. 

The magnetar population has increased very slowly 
since the discovery in 1979 of sources of r e- 
peated soft gamma-ray bursts (TMazets et al.l ll979aH bl: 
IGolenetskii et all 11984 ILaros et al.l 119861 I1987D . later 
associated with a new clas s of a s trophysical objects 
named SGRs (lAtteia et al.l 119871 : ILaros et al.l 119871 : 
iKouveliotou et al.lll987fr 7 and the confirmation, through 
bursting episod es, that AXPs we re part of the same 
group in 2002 ([Gavriil et al.l 120021) . New members are 
added in the group when (i) they are detected to emit 
multiple, soft short bursts and (ii) a spin period is 
found and a spindown rate is measured, which lead to 
magnetar-like -B-field estimates. During the 9 years 
of operation of the Compton Gamma Ray Observatory 
(CGRO; 1991-2000), we discovered only one new SGR 
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source, SGR 1627 - 41 (|Kouveliotou et all [1991 . In 
the first ~4 years of operation of the Swift satellite, no 
new source was discovered, althou gh several outbursts 
from known SGRs were recorded (iPalmer et al.1 120051 : 
llsrael et all 12005 lEsposito et all 120081 ) . The Fermi Ob- 
servatory was successfully launched on 2008, June 11 
and the F errm/Gamma-ray Burst Monitor (GBM) be- 
gan normal operations on July 14. During the first 
17 months of operation we recorded emission from four 
SGR sources, of which only one was a known magne- 
tar: SGR 1806 — 20. The other three detections were 
two brand new sources, SGR J0501 + 4516, discovered 
with Swift and extensively monitored with both Swift and 
GBM; SGRJ0418 + 5729, discovered with GBM, Swift 
and Konus-RF; and SGR J1550 — 5418, a source origi- 
nally known as AXP 1E1547.0-5408 or PSR J1550-5418 
(jCamilo et aLll2007l) . 

We present here the discovery of 
SGR J0418 + 5729 and its localization by triangula- 
tion in and the precise source localization with the 
Chandra X-ray Observatory and our infrared observa- 
tions in [J3] In section $4] we describe the properties of 
the GBM bursts and in ij5] we discuss the implications 
of our discovery. 

2. FERMI/GBM OBSERVATIONS AND LOCALIZATION BY 
TRIANGULATION 

The Fermi/ GBM has a wide field of view (8 sr) 
with broad-band energy coverage (8 keV — 40 MeV) 
and is, therefore, uniquely positioned t o detect tran- 
sient events (for a detailed description, see lMeegan et al.1 
I2009H . When GBM is triggered, three types of data 
are accumulated: CTIME Burst, CSPEC Burst, and 
Time Tagged Even t (TTE) data (jMeegan et al.l [20011 
IKaneko et all 120101) . The TTE data consist of time- 
tagged photon events for an accumulation time of ^330 s, 
starting ^30 s before the trigger time, with superior tem- 
poral resolution of 2 /us and fine spectral resolution of 128 
energy channels. 

GBM triggered on two SGR-like bursts on 2009 
June 5 at 20:40:48.883 UT and 21:01:35.059 UT 
(jvan der Horst et al.ll2009T ) . Their final on-ground calcu- 
lated locations, RA, Dec (J2000) = 70.0, +55.6 (4 h 40 m , 
+55°35') and 60.5, +55.4 (4 h 02 m , +55°22'), are shown 
in Figure [T] (top panels) . The positions are consistent 
at the ltr level with a common origin, and inconsistent 
at the 3cr confidence level with the known nearby SGR 
source, SGR J0501 + 4516, discovered with Swift in Au- 
gust 2008 (shown in the Figure Q] top panels). For both 
triggers, however, there is a systematic component to the 
localization uncertainty of 2 — 3° so that a reactivation 
of this known SGR could not initially be excluded, and 
indeed appeared the most likely origin for these events. 

The first GBM burst was seen also by the gamma- 
ray spectrometer, Konus-RF, onboard the CORONAS- 
PHOTON spacecraft and by the Swift/Burst Alert Tele- 
scope (BAT), which was triggered in its partially coded 
field-of-view. The second GBM burst was seen weakly 
in, but did not trigger, the BAT. Triangulation an- 
nuli of the GBM- Konus-RF and the GBM-BAT light 
curves for the first trigger are shown in Figure [T] 
(right upper panel); the GBM localizations are also dis- 
played. Subsequent ground analysis of the BAT data 
revealed a weak source at RA,Dec (J2000) = 64.606, 
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Fig. 1. — Top Left: GBM localizations (asterisks) for the two 
SGR J0418 + 5729 triggers on 2009 June 5. The contours indi- 
cate 1, 2, and 3<r statistical uncertainties. The position of the 
known SGR J0501 + 4516 is also shown. It is marginally con- 
sistent with the GBM events if one takes into account an addi- 
tional 2 — 3° systematic component to the localization uncertain- 
ties. Top Right: Triangulation annuli (90% confidence level) for 
the first SGR J0418 + 5729 trigger seen by GBM, Konus-RF and 
Swift/RAT on 2009 June 5. The asterisks and contours show the 
1<t (statistical uncertainty) localization of both GBM triggers (also 
shown in top left panel). A cross marks the position of the source, 
found in ground analysis using the Swift/HAT data (with a 4' 
uncertainty at 90% confidence level). Bottom: A 37" X 37" K B 
band image of the field of SGR J0418 + 5729 obtained with Palo- 
mar/WIRC. North is up and east to the left. Our CTiandra/HRC 
error circle with rad ius 0.35" is shown , as well as the original (3.6", 
Gogiis ct al. 2009b) and refined (1.9", Cummings et al. 2009) error 
circles obtained with Sm/i/XRT. 

+57.489 (4 h 18 m 25 s , +57°29'16") with an uncertainty 
of 4', which is consistent with the GBM localizations 
for both events (right upper panel in Figure [1]). The 
annuli clearly exclude the position of SGR J0501 + 
4516; the distance between the two sources is ~ 12°. 
Given these results we concluded that GBM detected 
SGR-like emission from a new source, which we named 
SGR J0418 + 5729 (|van der Horst et al.ll2009h . 

3. PRECISE SOURCE LOCATION WITH SWIFT, CHANDRA 
AND INFRARED OBSERVATIONS 

The Swift/X-R&y Telescope (XRT) observed 
SGRJ0418 + 5729 starting at 2009 July 8 at 20:52:35 
UT (when the source came out of Sun constraints for 
Swift) in photon counting mode for a total exposure 
time of 2.95 ks. A new X-ray source was found at 
RA = 04 h 18 m 33.70 s , Dec = +57°32'23.7" (J 2000) with 
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an error circle radius of 3.6". The XRT location is 
shown in Figure [1] (bottom panel) and is ~ 3.3' away 
from the initial Swift/BAT location. 

To obtain a precise location of the source, we observed 
the field containing SGRJ0418 + 5729 with the Chan- 
dra/High Resolution Camera (HRC) in imaging mode for 
23.8 ks on 2009 July 12. We constructed a binned (by a 
factor of two) image in the 0.5—7 keV band of the entire 
HRCT field and searched for point sources at 5cr above 
the background level. We discovered three previously un- 
catalogucd X-ray sources. We searched and detected co- 
herent pulsations at the spin period (9.08 s, IGogus et al.l 
l2009al) of SGRJ0418 + 5729 in the X-ray flux of the 
brightest X-ray source, which we id entified as the X-ra y 
counterpart of SGR J0418 + 5729 ([Woods et al.1 f2009h . 
The precise location (absolute positional uncertainty of 
0.35" at the 95% confidence level) of the source is shown 
in Figure Q] (bottom panel) and also given in Table [1] 
together with the coordinates of the other two sources. 

The Chandra position of SGR J0418 + 5729 was ob- 
served with the W ide field Infrared Camera (WIRC, 
I Wilson eT al. 2003) on the 5m Palom ar Hale telescope on 
2009 August T ifWachter et aDl2009h . WIRC has a field 
of view of 8.7' x 8.7' and a pixel scale of 0.2487"/pixel. 
We obtained 26 K s band images using two 13-position 
dither scripts (the second script spatially off-set from the 
first) - each position was a co-added exposure of twelve 
5 s images. Atmospheric conditions were very good with 
seeing < 1" and clear skies during observations. 

The individual frames were reduced using a suite of 
IRAF scripts and FORTRAN programs. These scripts 
first linearize and dark-subtract the images. A sky frame 
and flat field image are created from the list of input 
images, and subtracted from and divided into (respec- 
tively) each input image. At this stage, WIRC images 
still contain a significant bias that is not removed by the 
flat field. Comparison of 2MASS and WIRC photometric 
differences across the array shows that this flux bias has a 
level of ~10% and the pattern is roughly the same for all 
filters. Using these 2MASS-WIRC differences for many 
fields, one can create a flux bias correction image that 
can be applied to each of the "reduced" images. Finally, 
we astrometrically calibrated the images using 2MASS 
stars in the field. The images were then mosaicked to- 
gether and the mosaic was photometrically calibrated us- 
ing 2MASS stars. Vega magnitudes were computed using 
the IRAF phot routine with the zero points as found us- 
ing the 2MASS stars. The final image (see bottom panel 
of Figure [TJ has a 5a K s detection limit of 19.6 mag. 

The astrometric solution was derived based on 90 
2MASS source matches and carries a formal la error of 
0.1" for the transfer of the 2MASS reference frame to the 
WIRC image (in addition to the intrinsic 0.1" la uncer- 
tainty of the 2MASS reference system). The two addi- 
tional X-ray sources (see Table [T]) have unambiguous IR 
counterparts in our WIRC image and hence can be used 
to tie the X-ray astrometry to that of the IR 2MASS 
system. We find a small systematic offset between the 
two reference systems of 0.33" in RA based on those two 
sources and no systematic difference in Dec. The X-ray 
positions in TableQ]have been corrected for this shift and 
are thus registered to the 2MASS astrometric reference 
frame. 

Our K s band image overlaid with the X-ray error cir- 



TABLE 1 

Chandra/HRC coordinates OF SGR J0418 + 5729 and 

THE TWO X-RAY SOURCES WITH IR COUNTERPARTS 



Source 


Right Ascension 


Declination 


SGR.J0418 + 5729 


04 h 18 m 33 s .867 


+57°32'22.9l" 


(.12000) 


CXOU J041819. 0+573341 
CXOU J041812. 5+573154 


04 h 18 m 19 s .097 
04 h 18 m 12 B .578 


+57°33'41.75' 
+57°31'54.99" 


(J2000) 
(J2000) 



cles of IGogus et all (I2009bl ) the refined Swift/XRT po- 
sition bv [Cummin gs et al.l ()2009[ ) and our Chandra po- 
sition is shown in Figure [T] (bottom panel). No obvi- 
ous IR counterpart is detected inside the Chandra/HRC 
error circle. Two sources (labelled 1 and 2 in the 
Figure Q] bottom panel) with K s = 17.66 ± 0.04 and 
K s = 18.8 ±0.1, respecti vely, are detected within the re- 
fined Swift error circle of iCummings et all (|2009l ). How- 
ever, our Chandra/HRC position is sufficiently offset 
from this Swift/XRT position to exclude both of these 
sources. Possibly, a third, very faint source is seen at 
the southwestern edge of the Chandra/HRC error circle. 
Unfortunately, this source is at the detection limit with 
K s = 21.6±1.3 and cannot be reliably distinguished from 
a noise spike in the background. Hence our IR observa- 
tions fail to reveal a convincing counterpart candidate 
for SGR J0418 + 5729. 

4. SGR J0418 + 5729 BURST ANALYSIS 

We have searched the GBM continuous data files for 
untriggered bursts fro m SGR J0418 + 5729 using the al- 
gorithm described by iKaneko et al.l (|2010f ) starting two 
days before the two triggered events and ending six 
days after. Our search identified only three events from 
SGR J0418 + 5729, all detected on 2009 June 5: these in- 
clude one untriggered burst at 20:35:54.703 UT, and the 
two triggered events from the source. The untriggered 
event took place ~ 5 min before the first GBM trigger 
and is relatively weak and soft. It was detected only at 
energies < 50 keV, and its location is consistent at the 
la level with the SGR J0418 + 5729 location. We also 
checked the Swift/BAT data and although the source was 
in the BAT field-of-view, we did not see a rate increase 
at the event time. 

Additionally, a search through > 4000 IPN events with 
fluences > 7x 10~ 6 erg cm~ 2 and/or peak fluxes > 1 pho- 
ton cm~ 2 s _1 in the 25 — 150 keV energy range, going 
back to 1990, does not reveal a significant excess of bursts 
in the direction of SGR J0418 + 5729. We conclude that 
the source did not undergo any episode of intense ac- 
tivity during this time, although we cannot exclude the 
possibility of isolated, weak events similar to the ones 
reported in this Letter. 

We have performed detailed temporal and spectral 
analysis on the two triggered events using the TTE data 
type. The third (untriggered) event was only detected 
above background as one 256-ms bin in the continuous 
CTIME data. Since the CTIME data spectral resolution 
is relatively coarse (only 8 channels) , we did not perform 
a detailed spectral analysis for this event. 

The TTE data of the two triggered events were an- 
alyzed with the RMFIT (3.2rc2) spectral analysis soft- 
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Fig. 2.— TTE light curves of the two 
SGR J0418 + 5729 triggered events. The hatched areas indi- 
cate the time intervals used for spectral analysis. 



ware developed for the GBM data analysis 19 . We gener- 
ated Detector Response Matrices using GBMRSP vl.81. 
For both events we used in our analysis the three Nal- 
detectors with the smallest zenith angles (ranging from 
26° to 44°) to the source, i.e. Nals 3, 4 and 5. The 
Tgo and T50 event durations were estimated in RMFIT 
by constructing cumulative fluence plots over the en- 
ergy range 8 — 200 keV for the three detectors com- 
bined, and then determining the times during which 
90% and 50% of the b urst counts were accumulated 
(|Kouveliotou et ahl 119931 ). For the first trigger we find 
Tg {T 50 ) = 40 ± 7 ms (10±4ms), and for the second trig- 
ger T 90 (T 5 o) = 80 ± 6 ms (34±4 ms). 

Figure [5] shows the light curve in the detector with 
the smallest zenith angle for both events, with the in- 
tervals used for the spectral analysis indicated. We 
have fitted the time-integrated spectra with various func- 
tions: power law, cut-off power law, black body, and 
optically-thin thermal bremsstrahlung (OTTB). We find 
that OTTB provides the best fits in both cases (Table 
Ep similar to what has be en found for other SGR bursts 
(|G5gtis et alJ [l999. 2000J). The cut-off power law gives a 

19 R.S. Mallozzi, R.D. Preece, & M.S. Briggs, "RMFIT, A 
Lightcurve and Spectral Analysis Tool," ©2008 Robert D. Preece, 
University of Alabama in Huntsville, 2008 




Fig. 3. — Best fit OTTB count spectra of the two 
SGR J0418 + 5729 triggered events. 



better statistics value (Table[2]), which is not statistically 
significant given that this model has 1 additional free pa- 
rameter compared to OTTB. The best-fit count spectra 
are shown in Figure [3] From the Figure and the fit pa- 
rameters, it is clear that the spectrum becomes softer 
from the first to the second burst. 

To estimate the total energy output of the bursts, 
wc need to know the distance to the source. 
SGR J0418 + 5729 is located in the galactic plane and in 
the galactic anti-center direction. The biggest concentra- 
tion of stars in that direction is in the Pe rseus arm of ou r 
galaxy, at a distance of 1.95 ± 0.04 kpc (|Xu et al.H2006h . 
We provide here the SGR J0418 + 5729 burst energetics, 
assuming that the source is located in the Perseus arm, 
which we consider as an upper limit of the energetics. 
Adopting this distance implies energies (8—200 keV) of 
~ 4 x 10 37 and ~ 2 x 10 37 erg for the two bursts re- 
spectively, which is at the lowe r end of the distribution 
compared to other SGR bursts (iGogiis et al.lll999L I2000TI 
but at the high end for AXP ones (jGavriil et al.1 12004'). 
We note that an OTTB fit to the CTIME continuous 
data of the third (untriggered) event, implies an energy 
(8-200 keV) of - 8 x 10 36 erg, albeit with very low 
statistics. 

5. DISCUSSION 

During the first year of operations of the Fermi/GBM, 
we have detected bursts from four SGRs, of which two 
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TABLE 2 



Time 


OTTB 


Cut 


-off Power Law 




Energy Flux 




kT 


cstat/dof 


Index 




cstat / dof 


(8-200 keV) 


(UT) 


(kcV) 






(kcV) 




(10~ 6 erg/cm 2 /s) 


20:40:48.869 - 20:40:48.917 


33.46 ± 2.23 


296.76/361 


-0.51 ± 0.26 


34.72 ± 1.85 


294.85/360 


2.00 ± 0.08 


21:01:35.013 - 20:40:35.103 


19.71 ± 1.96 


335.15/361 


-0.66 ± 0.52 


21.39 ± 2.55 


334.85/360 


0.60 ± 0.04 



ones. SGR J0418 + 5729, in particular, was discovered 
with GBM and located by triangulation with Konus- 
RF and Swift. The source was subsequently confirmed 
as a magnetar candidate with th e Rossi X-ray Tim - 
ing Explorer (RXTE) observations (jGogiis et al.ll2009ah . 
which revealed a spin period of 9.0783±0.0001 sec in 
the persistent X-ray emission of the source. Although 
the source lies in the direction of the galactic anti-center 
and presumably at the nearby distance of ~ 2 kpc (as- 
suming it resides in the Perseus arm of our galaxy), 
it was not dete cted at any other w a velength (IR, op- 
tical and radio: iMignani et alJ 120091 : iRatti et al.l 120091 : 
lLorimer et al.l 120091 ) . Interestingly, this is the second 
source discovered in the same region of the sky (includ- 
ing SGR J0501 + 4516) in one year. The X-ray properties 
of the source's persistent emission will be described else- 
where (Woods et al., in preparation). 

The scarcity of magnetars contributes to the many 
open questions related to their nature, in particular the 
true number density and birth rate of these objects. The 
latest two new SGRs from roughly the same direction, 
if indeed at the relatively small distance of ^2 kpc, sug- 
gest that there are more such "dim" SGRs throughout 
our galaxy, undetectable unless they are relatively close 
to us. Indeed the trigger detection threshold of the two 
GBM triggers from SGR J0418 -I- 5729 indicates that if 
their origin was ~ 1.5 times further away, we would not 
have detected them. This raises the question: what is 
the population of such "dim" SGRs? A very rough esti- 
mate comparing 2 SGRs at ~2 kpc versus 4 at ~10 kpc, 
and assuming a uniform distribution within the galactic 
plane, gives ^10 — 15 times more "dim" sources active 
at a given time. To dominate the magnetar birth rate 
their active lifetimes should not be larger by more than 
this factor compared to those of "bright" SGRs. 

We have estimated the size of the parent population 
of SGRs using a technique that is commonly employed 
in the fields of biology and ecology to estimate animal 
populations (|Seberlll982l ). This ("Mark and Recapture") 
technique is based on capturing and marking a random 
sample of animals, and then returning them to the pop- 
ulation and allowing them to remix. When a new sample 
is captured later, the fraction of the recaptured animals 
that were already marked, and hence are in both sam- 
ples, can be used to es timate the p opulation size. We 
applied this technique ()Seberl [1982D to the SGR obser- 
vations assuming that the SGR population has (i) fixed 
membership, and (ii) is homogeneous in its bursting char- 
acteristics. There are several caveats associated with 
these assumptions (e.g., if some SGRs are quiescent for 
many decades, and then start bursting again, then the 
sample is biased towards a false "new" source), but we 
are using this estimate as a first order approximation. 



found by all instruments in the thirty years prior to GBM 
(note, that there was not always a complete sky coverage 
during this period and at times there were several large 
gaps, when no instrument was available to confirm SGR 
activity). The GBM observations have resampled the 
SGR population, finding 4 SGRs, 2 old and 2 ne w. The 
unbia sed form of the Lincoln-Petersen equation (|Seberl 
1982) estimates that the size of the SGR population is 
9 (+17.3/ - 1.6). The interpretation is that GBM is 
finding 50% old and 50% new SGRs, suggesting that the 
previously known sample is about half of the population 
observable by GBM. 

The discovery of SGR J0418 + 5729 adds a seventh 
confirmed member in the SGR subgroup of magnetars 
in the last 30 years. Together with the known AXPs, 
the total tally is ~ 15 sources, a very restricted member- 
ship club. Give n their small numbers, previous magnetar 
rate estimates (iKouveliotou et aD 11994}: Gacnsle r et alJ 
[pMlGTiTfe Hevlll2007HLeahv fc Ouvedll2009fl concluded 
that roughly 10% of neutron stars become magnetars. 
Our current rate estimates (based on the currently de- 
tectable sources) are c onsistent with the above, and with 
our earlier suggestion (|Kouveliotou et ail 1 1994) that our 
galaxy contains at any given time a few active magne- 
tar sources. However, if a dim, largely undetected as yet 
magnetar population exists, as the GBM detection of 
SGR J0418 + 5729 indicates, it might significantly con- 
tribute to and increase the magnetar birth rate. 
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